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We study the affinities between the shape of the bright soliton of the one-dimensional nonlinear 
Schroedinger equation and that of the disorder induced localization in the presence of a Gaussian 
random potential. With emphasis on the focusing nonlinearity, we consider the bound states of the 
nonlinear Schroedinger equation with a random potential; for the state exhibiting the highest degree 
of localization, we derive explicit expressions for the nonlinear eigenvalue and for the localization 
length by using perturbation theory and a variational approach following the methods of statistical 
mechanics of disordered systems. We numerically investigate the linear stability and "superlocal- 
izations". The profile of the disorder averaged Anderson localization is found to obey a nonlocal 
nonlinear Schroedinger equation. 



Introduction — Solitons and disorder induced An- 
derson states are two apparently unrelated forms of 
wave localization, the former being due to nonlinearity 
the latter to linear disorder [4]. However, on closer 
inspection, they look similar for various reasons: they 
are exponentially localized, they correspond to appropri- 
ately defined negative eigenvalues, they may be located in 
any position in space (which is homogeneous for solitary 
waves, and populated by a random potential for Ander- 
son states). Furthermore, various recent investigations 
deal with the theoretical, numerical and experimental 
analysis of localized states in the presence of disorder 
and nonlinearity [5l-fl9j , as specifically in optics 20 22], 



in Bose-Eistein condensation (BEC) |23H26j |. and more 
recently in for random lasers [27] . The nonlinear Ander- 
son localizations are expected to have a power (number of 
atoms for BEC, pump fiuence for random lasers or active 
cavities) dependent localization length, and eigenvalue, 
but also exist for a vanishing nonlinearity: in the low 
fiuence regime, they are Anderson localizations, but at 
high fiuence, it is expected that they are more related to 
solitons. Such a situation resembles other forms of linear 
localization, as the multidimensional "localized waves" 
|28], which are "dressed" by the nonlinearity 29|; a key 
difference with respect to localized waves is that Ander- 
son states are square- inte grab le, another feature in com- 
mon with bright solitons [30( . 

Many authors investigated the effect of nonlinearity on 
Anderson localization, as, e.g., 31 -35j], here we report on 
a theoretical and numerical analysis that allows to derive 
explicit formulae describing the nonlinear dressing of the 
fundamental Anderson states, and the way they become 
solitons as the nonlinear effects are dominant. We show 
that the disorder averaged profile of the nonlinear An- 
derson localization is given by the very same equation 
providing the soliton shape, augmented by a power and 
disorder dependent term. This equation defines a partic- 
ular highly nonlocal nonlinear response [3^ . and results 
in quantitative agreement with computations. In addi- 
tion, we numerically demonstrate that these states are 
stable with respect to small perturbations, and that this 



stability is driven by a novel kind of localization, which 
we address as "superlocalization" , resulting from the in- 
terplay of solitons and Anderson states. 
Outline — We review the nonlinear Schroedinger equa- 
tion with a Gaussian random potential; we describe the 
weak perturbation theory for small nonlinearity, when 
linear Anderson states are slightly perturbed; we consider 
the strong perturbation theory, i.e., the regime where 
the disordered potential is negligible and the only form 
of localization is the bright soliton; we compare the two 
limits with numerical simulations; we use a phase-space 
variational approach to derive results valid at any order 
of nonlinearity and quantitatively in agreement with the 
two mentioned limits and with numerical analysis; the 
stability is finally numerically demonstrated. 
The model We consider the one-dimensional 

Schroedinger equation with random potential V(x): 

i^t = -^ xx + V{x)iP- X \^\ 2 Tp = N[iP], (1) 

where x = 1 (x = — 1) corresponds to the focusing (de- 
focusing) case. V(x) is Gaussianly distributed such that 
(V(x)V(x')) =V ( ?6(x - x'). 

The linear states (x = 0) "0 = V ex P( — iEt) are given by 



Vxx + V(x)<p = Cip = Eip, 



(2) 



and we consider the lowest energy localized states with 
negative E; these are <p n {x) with (ipni^Pm) — S nm and 
S nm the Kronecker delta. We denote as E the lowest 
negative energy of the linear (x = 0) fundamental state 
ipo. We stress that ipo is the fundamental linear state 
with unitary norm, and in the following we will use states 
such that P = J | ip | 2 dx as the solutions of the nonlinear 
equation fl} ; P measures the strength of the nonlinearity 
as x — ±1 in the adopted scaling. 

The localization length is calculated by the inverse par- 
ticipation ratio 
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J ip^dx J (p^dx ^ 

For example, for an exponentially localized state ip^ = 
\/Pip e with ip e — exp (—2\x\/l)/y/l/2, we have I — I. Iq 



2 



is the linear localization length of the fundamental state 
y/Pipo. We recall that for the linear problem, the number 
of states per unit length is known and the mean value 
for the energy can be approximated by El — —Vq /3 



[4], l37| • Here El is the mean linear negative value of the 
energy levels of a Gaussian random potential, and will 
be used below as the appropriate limit for P — > 0. Wc 
also recall that in the linear case, the energy scales like 
the inverse squared localization length [H, as also found 
below when P — >• 0. 

The Lyapunov functional - The nonlinear Anderson 
states for a specific disorder realization V(x) are the so- 
lutions of 



]\[(ip) = -ip xx + V{x)ip - XV 3 = El P, 



(4) 



which are obtained numerically, as detailed below. The 
nonlinear states correspond to the extrema of the Lya- 
punov functional 

F = f m\ 2 + [V(x) - EM 2 - ^\ 4 }dx (5) 



with the Hamiltonian H = F + EP. 

As E = (ip,Af[<p\) /(if, <p) = (tp,J\f[ip\) /P, one has, for 
the solutions of (]!]), F = x{<P 2 , <p 2 )/2, that is 



F = ^ I <p 4 dx 



2 I ' 



1 2F 



I X P 2 ' 



(6) 



which show that a connection between the Lyapunov 
functional F and the localization length exists. 
Weak perturbation theory — For small P, standard per- 
turbation theory [l5| on y/~Pip (x) gives 



E = E -x^ + O(P 2 ), 



(7) 



where lo is the linear localization length in Eq.Q. For 
X = 1, E < decreases as the power is increased, while 
increases and eventually changes sign in the defocusing 
case (x = —1). Being ip^- 1 ' the standard first order cor- 
rection to the linear state ip , we find at order O(P) for 
the localization length: 
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(8) 

Eq.© predicts that I increases (decreases) with 
P in the defocusing (focusing) case; P = 
[^oY.nyo&n^l) 2 /{E n - E ))] gives the power 
level such that, when x = 1> I vanishes, and this is 
defined as the critical power for the transition to a soli- 
tonic regime, where the weak expansion is expected not 
to be valid. Pq depends on the linear eigenstates of the 
potential and comes from the lowest order perturbation 
expansion of the localization length. 



Summarizing, the weak expansion allows to affirm that 
two critical powers can be defined: in the defocusing 
case, there is a power P = \Eq\Iq at which the eigen- 
value changes sign, corresponding to a nonlinearity that 
destroys the Anderson states; in the focusing case there is 
a power P — Po at which the localization length vanishes, 
this is the nonlinearity level needed to the bound state 
for resembling a bright soliton (i.e., for the "solitoniza- 
tion" of the Anderson state). In the weak expansion these 
critical powers are dependent on the specific disorder re- 
alization and have a statistical distribution. In a later 
section, we report a variational approach that allows to 
derive closed expressions for the peak of these distribu- 
tions Pc, which depends only on Vb; we will limit to the 
focusing the defocusing one requires a separated 

treatment. 

Strong perturbation theory — For large P, we write the 
solution by a multiple scale expansion as ip = Pn(Px) 
and Eq.Q, at the highest order in P {xp = Px, E — 
(<p,Af[<p\) IP = P 2 E P ), reduces to 



V 



dxp 



(9) 



where Ep is the eigenvalue scaled by P 2 . For large P, the 
nonlinear Anderson states are asymptotically described 
by the solitary-wave solutions in a manner substantially 
independent of V(x). For x = 1, Eq.® is satisfied by 
the fundamental bright soliton Q and, correspondingly, 
we have cp = \/—2E/ cosh(^/~Ex). 

We stress that, in this expansion, E < as for the 
linear Anderson states, being E = E$ = — P 2 /16 and 
I = Is = 12/P, the subscript referring to the soliton 
trends. E$ is the "nonlinear eigenvalue" for the soliton, 
which is determined by P, while Is is the corresponding 
localization length in this strong perturbation expansion, 
i.e., when neglecting the linear potential V(x). Note also 
that this trend at high power is also expected for the 
higher order states <p n . 

Nonlinear dressing — We detail the transition from the 
linear Anderson states to the solitary wave. We limit 
to the focusing case X = 1 hereafter (the case x = — 1 
will be reported elsewhere). We numerically solve equa- 
tion Q, which is invariant with respect to the scaling 
x —} x/xq, V — ► Vxq and ip — > ip/xg, such that we can 
limit the size of system to x £ [— it, tt], over which period- 
ical boundary conditions are enforced. The prolongation 
of the linear states to the nonlinear case is not trivial. 
We start from a linear localized state (x = 0) and we 
prolong to x > by a Newton-Raphson algorithm; then, 
for each x, we rescale tp, by using the scaling properties 
of Eq.([T|), so that it corresponds to % = 1, and we cal- 
culate P, H, I and E. In figure [T^ we show the shape 
of the fundamental solution (lowest negative eigenvalue) 
for increasing power P. The inset shows the projection 
of the numerically retrieved nonlinear localization with 
the fundamental soliton sech profile: as P increases the 
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FIG. 1: (Color online) (a) Plot of the nonlinear Anderson 
states \ip\/max(\ip\) for different P (Vb = 4), two values of 
the powers are indicated corresponding to the blue and dashed 
lines; the inset shows the projection on the soliton profile; (b) 
nonlinear eigenvalue E versus P (Vb = 4) for several disorder 
realizations (cyan thin lines), compared with the strong (red 
thick line) and weak (green thick line, only shown for a single 
realization) expansions, and with the result from the annealed 
phase-space variational approach (black thick line). 
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FIG. 2: (Color online) (a) Localization length versus P 
(Vb = f): the results after various disorder realizations are 
shown (cyan thin lines) and compared with the weak (green 
line, only shown for a single realization with the correspond- 
ing Po indicated), with the strong expansion (red thick line), 
and with annealed phase-space approach (black thick line); 
(b) distribution of the critical power Po (200 disorder realiza- 
tions), the red vertical line is the analytical result for Pc, as 
given in the text (Vb = 2). 



shape of the disorder induced localization is progressively 
similar to the soliton. The trend of the eigenvalue E is 
shown in FigOJi, compared with the weak (low P) and 
strong expansions (high P), and with Eq. (fTT)) below: for 
low P we have a linear trend, Eq.([7]), while the trend 
follows the solitonic one E$ = — P 2 /16 for high P. 
In Figl^Ji we show the calculated localization length com- 
pared with the strong perturbation theory and, for a sin- 
gle realization, with the weak perturbation theory, with 
Po given by the intercept with the horizontal axis. When 
P increases, the localization length deviates from the lin- 
ear trend, and follows the bright soliton Is = 12/ P at 
high P for all the considered realizations. 
The phase- space/variational approach — In the weak ex- 
pansion Eq.©, valid as P —> 0, the power Pq for the 



transition to the soliton depends on the realization of 
the disorder and has a statistical distribution shown in 
FigHt". The histogram of Po m Fig(2}) is found not to 
substantially change for a number of realizations greater 
than 100. 

On the other hand, the strong expansion Eq.©, valid 
as P — > oo, completely neglects the random potential, 
and the result is independent of the strength of disorder. 
Here we introduce an approach based on the statistical 
mechanics of disordered systems .38] valid at any order 
in P. The first step is to define an appropriate measure 
based on the fact that the nonlinear Anderson state (HJ 
maximizes a weight in the space of all the possible ip. 
Following the fact the nonlinear bound states extremize 
the Lyapunov functional P, and that, for these states, F 
scales like l/l after Eq.©, we consider a Boltzmann like 
weight: Ti — exp(— L/l) with L determined in the follow- 
ing. Note that exp(— L/l) is the transmission of a slab of 
disordered material with length L and localization length 
I [3] ■ For a specific disorder realization, we introduce the 
measure 



pW\ = ^ ex P f-y 



with Z the "partition function" 



Z = 



exp 



dm, 



(10) 



(ii) 



such that J p[ip]d[ip] = 1. The inverse localization length 
is calculated as an average over the whole functional 
space of ijj [after Eq.©]: 



exp 



2LF 



d[# (12) 



In (fT2")) . the solution of (U]) is that providing the highest 
contribution to the weighted average among all the iji. 
Our aim is to find an equation for such a state after av- 
eraging over the disorder V(x); this averaging is denoted 
by an over-line: 



= —dhlog(Z) = —djJogZ. 



(13) 



In ([13]) w e used the so-called annealed average log(Z) = 
logZ,[39( whose validity is to be confirmed a posteriori. 
We find Z = J exp (-2LF eff /P 2 ) d[ip], being 



F, 



eff 



I^I 2 -tt(i + 



dx. (14) 



2 v- ■ p2 

This effective Lyapunov functional F e t t is extremized by 

2LV^ 



■ 4>i 



p 2 



\if>\ 2 ip = Ei/j 



(15) 



with the constraint P = J \ ip\ 2 dx, which gives E as a 
function P. Eg. (1151) generalizes the strong perturbation 



4 



limit Eq.©, retrieved for Vq — or P —> oo, to a finite 
potential Vq and P. 

Eq.([T5|) shows that the role of the disorder is to alter 
the nonlinear response, namely to increase the strength 
of the nonlinear coefficient, such that solitary waves are 
obtained at smaller power than in the ordered case. Con- 
versely, the linear localizations can be seen as the nonlin- 
ear Anderson states in the limit of vanishing power, that 
is a form of solitons only due to disorder. In Eq. (ji~5]l . P 
explicitly appears; this is due to the fact that the average 
over disorder introduces nonlocality (3(| in the model. In 
the defocusing case a result similar to (|T5"|) is found, with 
a nonlinear coefficient changing sign at high power, de- 
noting the absence of localization for large P, as we will 
report in future work. 

In the focusing case, by using the fundamental sech soli- 
ton of Eq (|15p . we find the corresponding eigenvalue, de- 
noted as Ec- 



E c 



p2 

16 



2LVq 

—p2~ 



(16) 
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FIG. 3: (Color online) Stability of the nonlinear Anderson 
states: (a) dashed, nonlinear state at power P = 7; full lines, 
two different superlocalizations arbitrarily vertically shifted 
(blu Q = 11, green Q = 14); (b) as in (a) for P = 23 (blu fi = 
59, green Q. — 61); (c) eigenvalues versus power P (Vq = 10); 
(d) evolution of the nonlinear states with P — 3.54 with 10% 
amplitude noise; (e) as in (d) with P = 13.2; the white line 
corresponds to the arbitrarily scaled initial profile tp (Vo = 2). 



with the localization length Iq = 3/ \J—Ec- Ec is the P 
and Vo dependent eigenvalue for a generic L, and lc is 
the corresponding localization length; note that accord- 
ing to this analysis a measurement of lc directly provides 
Ec- 

In the next step we determine L by imposing the cor- 
rect asymptotic value in the linear limit: as P — > 0, 

4/3 

it must be Ec — > El = —V a ' /3, which furnishes 
L = 2V i ^P/\fZ; conversely, in the large P limit one 
recovers the expected expression Ec — ► Es — — P 2 /16. 
Summarizing, we find the for the nonlinear eigenvalue 



Er 



p2 

16 



Pc 
P 



(17) 



with the only parameter Pc — 4V Q 2 ^ 3 Correspond- 
ingly, the localization length is 



la = 



12/ P 



(i+Pc/py 



(18) 



which also gives lg = 12/P for large P, and the weak 
limit Eq.© with l = 12/P C and P = P c . P C is the 
critical power for the transition from the Anderson local- 
izations to the solitons and is determined by the strength 
of disorder Vq- In the linear limit P — > 0, Eqs. (fT8| and 
(|17|) reproduces the known power-independent link be- 
tween the localization length and energy E = Ec — 
— 9/Zq JJ. We stress that l c is the localization length 
of the state that mostly contribute to Z. 
In Figures [T] and [2 we compare this theoretical approach 
with the numerical simulations at any value of P; results 
for various Vq are used to show that quantitative agree- 
ment is found in all of the considered cases. 



Stability and superlocalization — We consider the stabil- 
ity of the nonlinear localization: we calculate the eigen- 
values of the linearized problem following the Vakhitov 
and Kolokolov formulation [3^, 0, l4lj . We write tp = 
(ip + Sip) exp(— iEt), with Sip — [u(x) + iv(x)]exp(Ctt), 
where tp is a solution of Eq. @ and u and v are real val- 
ued. Eq.((3]) is linearized as 



fl u = LiLqu, 



(19) 



with the operators Lq = — d% — E + V(x) — fix) 2 and 
Li = -d'l - E + V(x) - 3(p(x) 2 . The stable (unstable) 
eigenvalues correspond to f2 2 > (f2 2 < 0). As it hap- 
pens for the standard solitons, the bound state profile is 
also an eigenvalue of p^|) . with 17 = due to the gauge 
invariance; conversely other neutral modes due to trans- 
lational invariance are lost due to the symmetry break- 
ing potential V(x). We numerically solve Eq.([T^l) and 
find that no unstable states are present for the consid- 
ered disorder realizations and values of Vo, demonstrat- 
ing that the nonlinear Anderson states are indeed stable 
with respect to linear perturbations. The interesting is- 
sue is that in regions far from the nonlinear bound states 
(where (p. = 0), Eq. ([T9| still admits non trivial solutions, 
corresponding to L\L = C 2 — [~d 2 + V(x) — E] 2 , such 
that the linear Anderson states correspond to f2 2 = 0. As 
P increases, the location of these states drifts towards the 
center of the nonlinear localization and this coupling re- 
sults into a power dependent O 2 (examples are given in 
Figj3^i,b,c). These can be taken as "superlocalizations" 
due to interplay between disorder and solitons. The sta- 
bility of the nonlinear Anderson states is also verified by 
their numerically calculated i-evolutions in the presence 
of a perturbation, as shown in Fig|3Ji,e. 



■5 



Conclusions — We reported on a theoretical approach 
on nonlinear Anderson localization demonstrating the 
strong connection between solitons and disorder induced 
localization. By a variational formulation we derived 
closed formulae for the fundamental state providing the 
trend of the nonlinear eigenvalue and the localization 
length at any power level in quantitative agreement 
with numerical simulations. Disorder averaged nonlin- 
ear Anderson localization is found to obey a nonlocal 
Schrocdingcr equation with a disorder dependent nonlin- 
earity. Such an equation, in the linear limit, reproduces 
the linear Anderson states. The reported approach can 
be extended to the multidimensional case and to other 
nonlincarities. 
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